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Abstract 

In a gas turbine engine, the turbine rotor blades are buffeted by the wakes of the vanes located upstream. There 
is a transient effect from the passing of wakes on the blade heat transfer. This transient effect has been computed for 
a representative rotor by introducing a wake upstream via an unsteady inlet flow boundary condition, or “gust” 
condition. Two cases of turbulent flow and laminar flow with Reynolds numbers of 385,000 and 385 respectively 
were considered. For the turbulent flow case a quasi-steady calculation was also performed. The variation in the 
unsteady heat transfer coefficient was found to be as high as 120 percent of the mean. For the turbulent flow case a 
quasi-steady calculation was also performed. The time mean of the unsteady heat transfer, the mean of the quasi- 
steady variations and the steady results agree reasonably well on all blade locations except for the turbulent results 
which differ near the leading edge. The quasi-steady heat transfer results do not agree with the instantaneous 
unsteady results, although the time-mean values are similar. 


Nomenclature 


C axial chord length 

h wall heat flux (Tw-Tt) 

k w fluid conductivity at the wall 

Nu Nusselt number = hC/k w 

p static pressure 

p 0 relative total inlet pressure 

Re Reynolds number based on the inlet velocity and axial chord 

S wetted surface distance measure from the geometric stagnation point 

negative on the pressure side and positive on the suction side 
Tt, Tw inlet total temperature and wall temperature 
U characteristic velocity, inlet axial velocity 

X axial location 

K> r reduced frequency = a>C/U 
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I. Introduction 


Turbine blades operate in an environment where the reduced frequency (roC/U) is outside of the small reduced 
frequency regime. The effect of the unsteadiness on the blade surface heat transfer and its significance is unclear. 
The effect of unsteady flow on heat transfer in other areas of turbine blades such as the blade tip and downstream of 
cooling holes are of interest. In the case of cooling holes, the time variation in discharge pressure forced by the 
passing of wakes causes the cooling stream to vary in time. The effect of this oscillation on the overall film-cooling 
effectiveness and the rate of heat transfer may be non-negligible and one which could be further investigated using 
unsteady CFD. The type of studies just described requires a three-dimensional unsteady calculation which our 
current work is planning to address. The work described herein is, however, limited to the surface heat transfer on a 
radial section of a blade subjected to wake induced buffeting, and is free from other complicating factors. 

The effect of flow unsteadiness on heat transfer has been studied under laminar flow conditions for a cylinder in 
an oscillating stream. The conclusion is that the effect of periodic unsteadiness, under the conditions of reduced 
frequencies that are much smaller than unity, can be neglected, and instead replaced by quasi-steady results. Quasi- 
steady conditions are defined by Lighthill as the state of the boundary layer that is appropriate to steady flow at the 
instantaneous value of the free stream conditions at the edge of the boundary layer. In effect, if the unsteady 
acceleration of oscillations is small this method can be utilized for simplification. Under the conditions of quasi- 
steady flow, described by Lighthill (ref. 1), the boundary conditions imposed on the boundary layer are not time 
variant. We have interpreted the quasi-steady condition as the condition for which one can specify the locations of 
the wake at the inlet and then solve the steady problem under those conditions. It should however be noted, that 
freezing the inlet boundary condition (i.e., keeping the wake from moving) results in the wake approaching the blade 
at a different relative angle. It will be seen later that this boundary condition still results in an average heat transfer 
that is almost identical to the time mean heat transfer but with the instantaneous value quite different than the 
unsteady value. 

The desire to obtain a complete unsteady flow and heat transfer solution in turbine blades is tempered by the 
enormity of the resources required. As the computational capabilities grow, obtaining unsteady solutions using 
turbulence closure models have become more realistic and indeed Unsteady Reynolds Averaged Navier-Stokes 
(URANS) solutions are becoming common. Two-dimensional solutions of unsteady rotor-stator interaction have 
been accomplished and such calculations have been used for clocking the blade rows (ref. 2) aimed at improving the 
efficiency of turbines. Unsteady RANS calculations including the rotor/stator interaction have also been performed 
to predict unsteady heat transfer on blades and vanes. Those would include the work of Rao et al. (ref. 3) who used a 
2 to 3 vane to stator blade count and a 2-D code and the work of Michelassi et al. (ref. 4) who rescaled the blade to 
maintain a one to one ratio of the vane-blade count and used a 3-D methodology to calculate the blade heat transfer. 
Abhari et al. (ref. 5) used the code UNSFLO which uses a transformed Euler equation to accommodate the vane to 
blade count ratio. The aerodynamics of passing of the wakes and the unsteady effect they have on low pressure 
turbine blades is dealt with in detail in a review paper by Flodson and Flowed (ref. 6). The main emphasis as related 
to low pressure blades is on the common problem of flow separation at condition of cruise speeds at high altitude. 
Under those conditions, the flow on the suction side of the rotor blade is laminar and prone to separation. The 
unsteady flow has the potential to be exploited to help energize the boundary layer and thus lessen the losses that 
would otherwise be produced by the separation. Fleidmann performed an analysis (ref. 7) and Fleidmann et al. 
performed an experiment (ref. 8) using unsteady wakes produced by a row of rotating rods passing over a 
downstream annular turbine blade cascade with showerhead film cooling. Unsteady cases were studied for a range 
of reduced frequencies as well as no-wake and quasi-steady wake cases. They found that increasing the frequency of 
rotation has a deteriorating effect on the effectiveness of film cooling. No heat transfer measurements were 
performed. In this paper we mainly address the issue of unsteady heat transfer on high pressure turbine blades. We 
will explore the suitability of quasi-steady calculations to replace the unsteady calculations as is at times done to 
reduce computational time. 
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II. Calculations 


Three-dimensional, URANS calculations were performed to predict heat transfer rate on a rotor blade. The flow 
was solved in an annular section assuming a slip boundary condition on boundaries possessing normals in the radial 
direction. Our simulation included a one to one and a two to one vane blade ratio with a gust type boundary 
condition at the inlet. The guide vane flow was calculated using the same code but in steady mode. The total 
pressure in the wake thus produced was fitted with a trigonometric function and placed as the boundary condition at 
the inlet of the rotor. The blade was rotating at 8400 lpm. There were 46 vanes and 76 blades (0.6 ratio). For the 
calculations presented here a one to one ratio was used. A one to two vane to blade ratio would be more accurate, 
but would double the computational cost and was not attempted. To study the effect of the one to one assumption, 
we doubled the frequency while holding the wake size the same to investigate the effect of frequency. The frequency 
of the wake thus produced is even higher than the actual configuration but was done to provide insight into the effect 
of frequency. 

The computer code used for this study was the Glenn-HT2000 (Steinthorsson (ref. 7)). This computer code was 
written using object oriented programming principles and the Fortran90 programming language. The numerical 
procedure uses a finite-volume discretization scheme that is second order accurate in time and space. The computer 
program uses the MPI (Message Passing Interface) parallel processing. The cases were run on 20 processors of a 98 
processor Xeon Linux Cluster. The turbulence model used for the calculations was the Low Reynolds number k-co 
model of Wilcox integrating to the walls. The grid used was quite fine adjacent to the blade (y + ~l), in the boundary 
layer and in the free stream. The turbulence model is able to produce an effect similar to the transition from laminar 
to turbulent flow. In practice however, the transition is not guaranteed to be in the appropriate location. In fact, it is 
often triggered very near the leading edge which is what occurred in our computations making them fully turbulent 
except near the leading edge. The following was investigated: 

1. Unsteady effect of the passing wake on the rotor blade row downstream for a laminar and two turbulent cases at 

0.44 and 0.88 reduced frequency. Unsteady and time average. 

2. Steady calculation particular to the average inlet condition of the unsteady case. 

3. Quasi-steady flow calculations and determination of the ensemble average of the blade heat transfer. 

4. Laminar flow calculation of heat transfer over the blade. 

The quasi-steady case was performed as a series of 32 steady calculations with the wake profile specified at equally- 
spaced intervals across the inlet. This was done to determine if the motion of the wake passing is important relative 
to the simple presence of the wake. In effect, the quasi-steady case is equivalent to the limiting case of zero reduced 
frequency. 

The unsteady solution process was started from a converged steady initial solution and continued for several 
cycles until a transient periodic solution to the blade heat transfer was achieved. The final time step was chosen by 
progressively reducing the time step until there was little change in the unsteady results. The convergence efficiency 
was also considered and it was then decided to cover each period by 640 time steps. 

The particular geometry considered for this case was the blade geometry and flow conditions of the GE-E 3 
(ref. 10). A simplified case and a full three-dimensional geometry were considered. The three-dimensional 
calculations are still in progress and will be presented elsewhere. The simplified geometry was a thin ‘sliver’ of the 
blade for which slip boundary conditions in the upper and lower radial boundaries were specified. There were four 
grid intervals in the radial direction and the blade-to-blade grid was constructed as shown in figure 1 . The grid was 
refined to resolve the passing of the wake. A coarse grid does not support the wake and the wake dissipates quickly. 
The simplified case was set to rotate about the rotation axis and was also subjected to a gust type boundary condition 
upstream. This case was used to investigate the effects described above. The blade row pressure ratio was 0.44 and 
rotation rate was 8400 rpm. The case was scaled to give a Reynolds number based on the axial chord and inlet 
velocity of 385,000 for the turbulent flow cases. The Reynolds number for the laminar flow case was 385. The inlet 
gust boundary condition was specified such that it has a 15 percent total pressure deficit. It has a background 
turbulence intensity of 2 percent and reaches an intensity of 3.5 percent in the wake. These values were arrived at 
after the flow through the vane row upstream was calculated independently. 
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Figure 1 . — Blade profile and the grid in a radial plane. 


Results 



Figure 2, shows the time variation in the local 
static pressure on the blade surface. This variation in 
the surface pressure would have repercussions if 
film cooling were performed on the blade. In fact, it 
would have resulted in a pulsation of the film- 
cooling flow. In the figure, in addition to the time 
average of the pressure we are showing the pressure 
distribution resulting from a steady calculation. The 
difference in the steady solution and the time 
dependent solutions with the wake is obvious. The 
difference is attributable to the difference in the 
incidence angle which would be less steep, on the 
average, due to the presence of the wake. In other 
words, because we hold the absolute angle as fixed, 
the reduced velocity in the wake produces a reduced 
relative angle of attack. This is an important point to 
consider when running a CFD code for a rotor in 
steady mode. 

Figure 3 shows the unsteady heat transfer under 
laminar flow conditions. The time averaged values 
and steady results are also shown in that figure. 

Those differences are quite small. As the flow is 
purely laminar, the unsteady variation in heat 

transfer has to be due to flow acceleration. The correlation coefficient between the unsteady components of shear 
and Nusselt number was computed and was found to have a value near unity. The correlation was computed taking 
into account all the points on the blade surface over a complete period. The interpretation is that the heat transfer 
increase/reduction is due to flow acceleration/deceleration under low Reynolds number conditions. The 


Figure 2. — Relative 
distance from the 
0)r= 0.44. 


pressure distribution 
leading-edge. Re = 


versus axial 
385,000 and 
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Figure 3. — Nusselt Number versus wetted 
distance for laminar flow. Re = 385, and co r = 
0.44. 


Figure 4. — Nusselt Number versus wetted 
distance for turbulent flow calculations. Re = 
385,000 and w r = 0.44. 


good correlation suggests that there is not a large lag 
between these two effects in laminar flow situations. 

Figure 4 is a plot of Nusselt number for turbulent 
flow over the rotor blade, the difference between the 
average of unsteady calculations over a period and 
steady calculations does not appear to be large 
except near the leading edge. This may be due to the 
differences in the incidence angle or the unsteady 
transition effect or both. An analysis of the unsteady 
components of shear and Nusselt number yielded a 
small correlation coefficient. This assertion can be 
verified by noting that the mean heat transfer and 
steady value are matched near the leading edge for 
the laminar flow case. The effect was explained by 
Dullenkopf and Mayle (ref. 11) who also proposed a 
transition model which help predict the time average 
heat transfer on the blade. This model was used in a 
calculation of heat transfer over a rotor blade. 

(Ameri and Arnone (ref. 12)). The real effect is not 
fully apparent here since the model utilized in this 
work does not include a physically based bypass 
transition model. In many instances, such as 
instances in which there is film cooling near the 
leading edge, or large enough adverse pressure 

gradient, transition does occur early. The relative magnitude of the envelope obtained for the laminar and the 
turbulent flow calculations in figures 3 and 4 is shown in figure 5. The figure shows that the unsteady variation 
around the mean can be as large as 1 10 percent of the mean. It is interesting that the unsteady variation experiences 
a precipitous drop once it gets to the uncovered portion of the blade. 

In figure 6 the vector plots correspond to the unsteadiness (instantaneous — time mean) at two different phases of 
the periodic unsteady wake passing for the turbulent flow computations. The shaded background is the entropy rise. 
The lighter areas correspond to the wake locations. 



Figure 5. — Percent Variation in Nusselt 
NuMin)/NuAvg versus axial distance from 
edge. 


NO. (NUwax- 
the leading- 
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Figure 6. — View of unsteady component of velocity distribution at two different snapshots in time 
Re = 385,000 and co r = 0.44. 


Figure 7 shows the unsteady heat transfer for the 
same flow as in figure 4 except that the wakes occur 
twice as often. The motivation for this exercise was 
explained earlier. The figure shows that the size of 
the unsteady heat transfer envelope has been 
reduced considerably compared to figure 4. As the 
case in figure 4 was run with a one to one ratio of 
vanes to blades (itself too fast by a factor of 1.6), it 
may therefore be concluded that running with the 
proper ratio of vanes to blades would likely result in 
even a larger unsteady heat transfer envelope. 

Quasi-Steady Results 

There have been instances in which workers 
have used the quasi-steady conditions to simulate 
unsteady conditions. As Lighthill (ref. 1) points out 
this condition would be a good substitute for the 
unsteady case as long as the reduced frequency is 
small. In order to obtain quasi-steady solutions the 
wake was moved on a number of locations, of equal 
spacing, along the inlet boundary and corresponding steady solutions was obtained for those locations. The solution 
envelops and the average values are also shown in figure 8. It is observed that the envelope of the unsteady heat 
transfer (for the lower frequency), in figure 4, is much wider than that of the quasi-steady solution on the suction 
side. The opposite is observed on the pressure side. The average however, is very close to the average heat transfer 
obtained for unsteady flow. The difference points to the fact that the quasi-steady solution cannot replace the 
unsteady solution. But, it is interesting that the cases so far in this work show that either the quasi-steady mean or 
the steady value can be substituted for the time-mean value. This is especially true on the suction side. On the 
pressure side, some differences do exist but, the differences may be too small. The one location where the 
differences appear to exist and could be significant is near the leading edge of the blade. Both the quasi-steady 
average and steady values are significantly different than average. Figure 9 is the quasi-steady counterpart to 
figure 6. 



S 

Figure 7. — Higher frequency, Nusselt number versus 
surface wetted distance. Re = 385,000 and co r = 0.88. 
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s 

Figure 8. — Quasi-Steady flow, Nusselt number 
versus surface wetted distance. 



Figure 9. — View of deviation of the velocity vector 
from the mean for a quasi-steady frame, 
Re = 385,000. 


NASA/TM- 


-2007-214908 


7 



Figure 10. — Nusselt number variation at a point on the 
suction side of the blade located at S = 1, 
downstream of the leading edge for unsteady and 
quasi-steady computations. 

Figure 10 shows the phase variation of Nusselt number for the two blade count ratios considered and the quasi- 
steady solution at those phases for a location one chord length downstream of the leading edge on the suction side of 
the blade. The figure shows that the quasi-steady values should not be used in place of unsteady values. 


Conclusions 

In this paper the effect of unsteady wake in a simplified model of the rotor/stator interaction was considered. The 
wake was modeled after the vane flow was independently simulated. The results are limited to larger reduced 
frequencies as compared to those that occurred in the real situation because only one blade periodic case was 
simulated. Two frequencies were considered. It was found that: 

1. The steady heat transfer results are good representation of the mean heat transfer for a clean blade (no film 
cooling and 2-D flow) for all the cases considered. 

2. An average computed from a quasi-steady flow is also a good substitute for time average Nusselt number. 

3. The quasi-steady value is not a good representation of the heat transfer on the blade for the unsteady values and 
in fact the envelopes are substantially different. 

4. An increase in “reduced frequency” of wake passing from 0.44 to 0.88 (doubling) resulted in a substantial 
shrinking of the unsteady heat transfer envelope. 
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